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A  numerical  investigation  of  head-on  collisions  of  solitary  waves  in  water  shows  that  they  reform 
to  be  solitary  waves  almost  as  before  the  collision,  but  not  quite.  A  little  energy  is  deposited  in  a 
dispersive  wave  train  generated  during  the  collision.  Hence,  solitary  waves  in  water  are  not  exactly 
solitons.  The  calculations  further  indicate  a  transition  where  the  water  at  the  crest  attempts  to  fly 
off  the  surface  when  an  identical  pair  of  incoming  waves  have  an  amplitude /depth  ratio  exceeding 
approximately  0.45. 
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SOLITARY  WAVES  IN  WATER  COLLIDING  HEAD-ON 


When  solitary  waves  in  rater  collide  head-on,  a  brief  but  nonlinear 
interaction  occurs,  after  which  the  solitary  waves  refora.  If  exactly 
unchanged  in  fora,  the  use  of  the  tern  "soli ton, "  coined  by  Zabusky  and 
Kruakal*,  is  appropriate.  Soli tons  are  solutions  of  many  model  equations 
that  mathematically  describe  slightly  nonlinear  waves  in  many  media,  the 
first  equation  shown  to  possess  soli ton  properties  is  the  Korteweg-deVries 
Equation2,  originally  derived  to  describe  the  evolution  of  rater  raves. 

Like  Ref.  2,  this  report  deals  with  water  raves.  It  addresses  the  question 
of  how  closely  solitary  raves  following  a  head-on  collision  resemble 
soli tons. 

Canfield  and  Street3  and  Maxworthy4  studied  these  collisions 
experimentally,  by  reflecting  a  solitary  rave  from  a  rigid  wall  or  by 
running  two  solitary  raves  into  each  other,  they  find  that  the  maximum 
amplitude  during  the  collision  is  larger  than  the  sun  of  the  initial 
amplitudes,  by  an  amount  that  increases  with  increasing  amplitude. 

Maxworthy  reports  that  a  new  phenomenon  occurs  when  the  amplitudes  of  the 
colliding  raves  become  sufficiently  large,  the  rater  near  the  crest 
attempts  to  lift  off  the  body  of  water  and  subsequently  collapses  back  into 
the  rest  of  the  rater.4  For  a  pair  of  equal-amplitude  solitary  waves 
meeting  inside  the  rave  tank,  the  transition  occurs  when  the  solitary  raves 
start  with  amplitudes  near  0.5.  When  a  rave  reflects  from  a  wall,  the 
effect  is  deferred  to  larger  amplitudes,  probably  because  of  cohesive 
forces  between  the  rater  and  the  wall. 

the  theoretical  models  of  colliding  solitary  raves  take  the  rater  to  be 

inviscid  and  incompressible,  and  the  flows  to  be  irrotational.  Chan  and 

Street3  applied  a  modified  marker-in-cell  numerical  method  and  report  "In  a 
Msasistlpt  approved  Sspfmbsr  14, 1983. 


key  test  of  the  method,  a  wave  was  reflected  from  the  wall  and  returned  to 
Its  initial  position  for  the  case  HQ/d  »  0.2.  [RQ  is  the  initial 
solitary  wave's  elevation  above  still  water  depth  d. J  The  reflected  wave 
had  exactly  the  same  surface  profile  as  the  incident  wave." 

Analysis  has  been  applied  to  this  problem,  most  recently  to  third  order 
in  ff0/d  by  Su  and  Mirie.8  They  find  that  to  this  order  the  waves 
asymptote  to  solitary  waves  identical  in  form  to  the  initial  waves,  though 
sliqhtly  retarded  in  location.  They  also  predict,  however,  the  existence  of 
a  dispersive  wave  train  generated  during  the  collision.  The  dispersive  wave 
train  has  some  energy  (not  necessarily  at  third  order)  that  has  to  come  at 
the  expense  of  the  solitary  waves.  Conseguently,  the  existence  of  a 
dispersive  wave  train  guarantees  that  the  outgoing  solitary  waves  are  less 
energetic  than  the  incoming  solitary  waves.  In  short,  the  solitary  waves 
are  not  exactly  soli tons. 

Fenton  and  Rienecker7  recently  reported  the  application  of  an  accurate 
numerical  method8  to  colliding  solitary  waves.  They  report  substantial 
changes  in  solitary  wave  amplitudes  after  collision,  and  cast 
doubts  on  the  analysis  of  Su  and  Mirie.  They  do  not  find  evidence  of  a 
dispersive  wave  train  throuoh  the  admittedly  brief  duration  of  their 
numerical  experiments. 

This  report  describes  the  application  of  the  "unified  waves  model”9  to 
head-on  solitary  wavs  collisions.  This  model  uses  finite  difference 
technigues  for  propagation  in  one  horisontal  dimension,  it  is  capable  of 
treating  both  irrotational  and  rotational  waves  in  channels  of  gradually 
varying  depth  and  breadth,  with  ambient  steady  currents  present,  and  with 
forced  or  free  waees.  Model  results  have  been  compared  with  results  from 
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laboratory  experiments  and  with  those  from  other  numerical  and  analytical 
theories.  10  date,  the  comparisons  with  theory  show  the  model  to  be  very 
accurate;  some  features  seen  in  experiments  are  seen  in  the  unified  waves 
model,  but  not  in  other  models.  Features  of  the  model  are  summarized 
below. 

1.  The  model  uses  exact  prognostic  equations  in  conservation  form.10  These 
appear  not  to  have  been  used  in  wave  modeling  before  now.  An  integral  of 
the  momentum  equation  is  an  extension  of  Kelvin's  circulation  theorem. 

2.  Higher  order  expansions  than  used  in  other  numerical  models  connect  the 
velocity  variables  that  appear  in  the  governing  equations.  This  allows  the 
model  a)  to  incorporate  long  wave  theory  exactly,  b)  to  include  both  shallow 
water  and  essentially  deep  water  waves  in  the  same  model,  and  c)  to 
represent  fairly  long  nonlinear  waves  to  one  order  better  than  Boussinesq  — 
hence  the  name  "unified  waves  model." 

3.  The  model  employs  a  numerical  method,  pure  leapfrog,  that  gives  no 
unwanted  numerical  diffusion.  The  time-stepping  procedures  are  simple 
enough  to  analyze  in  some  detail  and  to  implement  efficiently  on  vector 
computers . 

4.  The  model  can  take  a  time  step  equal  to  a  space  step  (in  nondimens iona 1 
units  in  which  the  linear  long  wave  speed  is  unity).  This  permits  efficient 
machine  computations ,  unlike  methods  developed  for  the  Korteweg-devries 
Equation.  Moreover,  this  procedure  removes  any  spurious  numerical 
dispersion  at  the  lowest  order. 

5.  Finally,  the  diagnostic  equations  are  cast  in  a  form  such  that  only 
tridiagonal  matrix  equations  need  to  be  solved.  A  very  fast,  fully 
vectorised  algorithm  is  then  used  to  invert  the  matrices. 
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The  first  three  columns  of  Table  1  list  properties  of  the  solitary 
waves  used  here*  These  are  obtained  by  running  an  initial  disturbance  that 
resembles  a  solitary  wave  through  the  computational  channel,  allowing  the 
solitary  wave  to  outrun  the  rest  of  the  disturbance,  and  extracting  the 
solitary  wave  from  the  rest  to  use  as  initial  conditions  for  runs  involving 
collisions.  The  solitary  wave  properties  given  in  Table  1  are  very  close  to 
those  of  exact  solitary  waves,  as  derived  by  Longuet-Higgins  and  Fenton.1'1' 
The  accuracy  is  about  that  of  a  fourth  order  expansion, ^  though  the  formal 
accuracy  is  limited  to  second  order  in  the  treatment  of  the  diagnostic 
equations.  Perhaps  the  unwarranted  accuracies  are  due  to  retaining  exact 
prognosis  thoughoutj  perhaps  they  are  fortuitous. 

Figure  1  provides  a  global  view  of  the  collision  of  two  solitary  waves, 
each  having  an  initial  amplitude  of  0.386.  The  solitary  waves  appear  to 
survive  the  collision  intact,  weak  additional  disturbances  are  visible  in 
the  region  between  the  solitary  waves  after  the  collision,  however. 

Figure  2  employs  the  same  data  as  Fig.  1,  but  increases  amplitudes  10- 
fold,  and  clips  the  profiles  at  0.02.  Figure  2  brings  out  some  details  of 
the  additional  disturbance.  Clearly,  it  is  a  dispersive  wave  train 
generated  by  the  collision,  consistent  with  the  prediction  of  Su  and  Mirie . 

Figure  3  plots  the  amplitude  and  phase  of  the  wavs  crest  that  moves 
left  to  right  in  Figure  1.  The  behavior  is  characteristic  of  all  runs. 

After  a  brief  transient  the  crest  reaches  a  constant  value,  merges  with  the 
other  solitary  wave  (at  t  s  100),  where  the  amplitude  more  than  doubles, 
asperates,  reaches  a  minimum  amplitude  (just  before  t  ■  120),  and  slowly 
reforms  to  a  solitary  wave  slightly  less  high  than  before  the  colliion. 
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The  phase  curve  increases  at  first,  indicating  an  accelerating  crest 
until  the  two  crests  marge  at  x  -  200.  The  crest  renains  in  this  location 
for  a  while,  and,  upon  the  reemergence  of  two  crests,  accelerates  until 
t  3  140,  after  which  a  crest  decelerates.  Because  the  asymptotic  outgoing 
waves  travel  a  little  slower  than  the  incoming  ones,  the  phase  lag  will 
increase  as  time  goes  on.  By  extrapolating  the  phase  line  shorn  in  Fig. 3 
from  the  region  between  t  ■  210  to  t  »  268  back  to  the  collision  time 
t  -  56.675  one  obtains  an  estimate  of  the  phase  lag.  This  turns  out  to  be 
&  -  0.478. 

Numerical  tests  were  performed  to  determine:  1)  Can  truncation  error 
be  generating  the  dispe naive  wave  train,  and  so  be  responsible  for  the 
slightly  non -soli ton  behavior?  2)  Is  the  particular  second  order 
formulation  chosen  accidently  giving  results  that  are  not  common  to  other 
formulations  accurate  to  the  same  order?  3)  Can  coding  errors  be 
introducing  artificial  diffusive  or  anti -diffusive  effects?  4)  Are  the 
quantities  that  should  be  conserved  in  the  calculations  actually  conserved? 

Runs  at  different  resolutions  Ax  were  performed  to  answer  question  ( 1 ) . 
Specifically,  for  head-on  collisions  between  solitary  waves  having  initial 
amplitudes  of  0.386  were  run  with  three  resolutions,  &c  ■  1/4,  Ax  ■  1/8,  and 
Ax  «■  1/16,  each  with  At  -  Ax.  The  results  differ  slightly,  particularly  in 
the  phases  of  the  crests  and  troughs  of  the  dispersive  wave  train.  Yet  if 
truncation  error  is  responsible  for  the  dispersive  wave  train,  one  expects 
amplitudes  to  be  proportional  to  (  Ax)  z,  or  to  vary  by  a  factor  the  order  of 
16  between  the  Ax  ■  1/4  and  Ax  ■  1/16  runs.  The  observed  amplitude 
variation  is  only  a  few  percent.  Consequently,  we  can  rule  out  truncation 
error  as  being  responsible  for  the  existence  of  the  dispersive  wave  train. 
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ft*  unified  waves  nodal  is  flexible  in  evaluating  high-order  linear 
dispersive  terns.  Tor  problem  for  which  nonlinear  and  dispersive  effects 


are  in  approximate  balance,  such  as  the  problem  at  hand,  the  theory  is 
probably  accurate  only  to  second  order,  no  matter  how  accurately  the 
dispersive  tens  have  been  taken  into  account.  Runs  with  fourth  order 
dispersivt  terms  included  agree  very  well  with  runs  that  drop  dispersive 
terns  beyond  second  order.  Consequently,  question  (2)  above  is  answered,  in 
that  at  least  two  partially  independent  second-order  formulations  give  the 
sane  results . 

The  formulation  and  numerical  methods  introduce  no  numerical  diffusion 
or  antidiffusion.  TO  check  whether  any  coding  errors  introduce  diffusion,  a 
run  that  involved  time  reversal  was  performed.  For  this  run  Ax  -  At  ■  1/8, 
and  the  initial  solitary  wave  started  near  x  -  131,  bounced  off  a  wall  at 
x  -  200,  and  returned  to  x  -  152  over  fit  -  100  (800  time  steps).  Time  was 
then  reversed  and  the  run  continued  for  another  800  time  steps  (fit  «  -100). 
Rssults  are  compared  at  time  steps  50  and  1550.  With  a  perfect  computer, 
the  fields  at  these  time  levels  should  be  identical.  The  elevation  and 
surface  velocity  at  fixed  grid  points  differ  slightly,  but  by  no  more  than 
0.000024  in  elevation  and  0.000061  in  surface  velocity.  These  deviations 
are  consistent  with  roundoff  error.  Consequently,  diffusive  coding  errors 
axe  absent  (or  result  in  unde tec tab ly  tiny  effects) . 

The  final  question,  that  of  proper  conservation,  is  a  bit  complex.  The 
formulation  and  numerical  implementation  are  designed  to  conserve  two  linear 
field  quantities i 

/  ndx  ■  Constant  +  Hass  injection  from  ends  (1) 
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/  (u  +  n'v  )dx  ■  Constant  +  Momentum  injection  from  ends  (2) 

s  s 

In  (1-2)  n  is  the  elevation  above  still  water  level,  u_  and  v  are  the 

s  s 

components  of  horizontal  and  vertical  velocities  of  the  fluid  at  the  free 
surface,  and  n'  is  the  surface  slope.  The  integrals  span  the  computational 
channel,  left  to  riqht.  ft?.  (1)  is  the  obvious  conservation  of  mass; 
eo.  (2)  is  an  extension  of  Kelvin's  circulation  thorem1®.  The  numerical 
methods  are  designed  to  satisfy  (1)  and  (2)  exactly,  and  constant  monitorino 
of  both  sides  of  (1)  and  (2)  show  that  the  methods  work  perfectly  to  the 
four  significant  figures  monitored.  Hie  energy,  however,  which  the  physics 
demands  be  conserved,  is  not  conserved  exactly  by  the  unified  waves  model, 
the  calculations  add  some  energy  to  the  computational  domain  during  a 
solitary  wave  collision,  and  remove  most,  though  not  all,  of  the  energy  as 
time  goes  on.  the  net  energy  addition  leaves  some  uncertainty  about  the 
change  of  solitary  wave  amplitudes.  If  one  takes  away  the  artificially 
introduced  energy  from  the  emerging  solitary  waves,  they  would  have  slightly 
smaller  amplitudes.  Using  interpolations  from  table  5  of  the  paper  by 
Longuet-Riggins  and  Fenton11  one  can  relate  solitary  wave  amplitude  and 
total  energies.  I  regard  the  failure  of  the  unified  waves  model  to  exactly 
conserve  energy  as  introducing  an  uncertainty  in  the  amplitude  of  the 
emerging  solitary  waves,  but  not  in  the  general  conclusion,  to  be  drawn, 
that  the  colliding  solitary  waves  are  not  guite  solitons. 

Some  Quantitative  results  of  the  collisions  of  two  solitary  waves  with 
themselves  and  with  each  other  are  given  Table  1.  One  property  of  solitary 
wave  collisions  not  reported  in  other  theoretical  works  is  the 
maxima  vertical  acceleration  at  the  surface.  This  exceeds  that  of  gravity 
in  the  calculations  of  the  0.4726  collision.  This  would  imply  negative 
pressures  inside  the  fluid  near  the  surface,  and  the  fluid  there  would 


attempt  to  "fly  off.  Mora  precisely,  it  would  bo  loft  behind  by  the  faster 
falling  fluid  below  it.  Thus,  new  phenomena  not  deacribable  by  the  unified 
wawee  aodel  are  starting  when  initial  solitary  wave  amplitudes  are  less  <■»>».» 
0.4726.  Presumably,  these  new  phenomena  are  those  found  experimentally  by 
Hexworthy.  **  A  special  set  of  calculations  was  performed  to  identify  the 


point  at  which  the  maximum  vertical  deceleration  matches  g.  The  initial 
conditions  involved  solitary  waves  outrunning  a  residual  disturbance. 
Calculations  with  he  »  At  -  1/16  were  run  in  a  computational  domain  of  3200 
grid  points  through  the  time  that  the  crest  reached  its  maximum.  The 

calculations  indicate  that  the  transition  should  occur  at  initial  amplitudes 
of  0.45. 

Bis  other  results  of  htble  1  agree  well  with  the  results  of  at  least 
one  previous  investigation.  The  behavior  during  the  collision,  as  measured 
by  the  comparison  of  runup  with  FWnton  and  Rienecker,  and  with  experiments, 
i»  contorting,  because  one  cannot  guarantee  good  performance  from  the 


unified 


Model  for  the 


amplitudes  encountered  during  the 


collision.  The  substantial  decrease  in  wave  amplitude  at  the  early  times 
eaantne*  by  Kenton  and  Rienecker  is  confirmed. 

Bis  eventual  relaxation  to  solitary  waves  only  slightly  lees  energetic 
then  ini  tally,  SM,  especially,  the  strong  dependence  of  da  on  amplitude 
(about  gnartic)  are  consistent  with  the  null  result  te  -  o  given  from  the 
third  order  theory  of  flu  and  Riria.  These  small  changes  must  be  real, 
hiaNnfar#  as  Figurs  2  and  energy  conservation  neks  clear.  Thus,  the  term 


in  water,  for  they 


solitary 


tbs  'sms  attar  colliding  haad-oa,  though  almost. 


1.  All  computations 


All  entries  are  in  units  g  ■  HQ  ■ 
that  generate  the  data  here  take  £x  ■  At  ■  1/8.  The  first  three 
columns  list  solitary  wave  properties  that  have  evolved  to  a 
constant  value  prior  to  collision.  The  figures  in  brackets  use 
expansions  or  table  interpoltions  from  Longuet-Higgins  and 

m 

Fenton.11  Solitary  wave  mass  is  /  ndx. 

“«• 

Columns  4  and  5  are  measured  at  the  time  of  maximum 
elevation.  "Runup"  is  the  elevation  above  still  water  level  and 
•Accel,"  is  the  magnitude  of  the  vertical  acceleration  of  the 
fluid  particle  at  the  crest.  Here,  and  in  the  other  entries,  the 
numbers  in  ordinary  parenthases  are  from  the  third  order 
analytical  results  of  Su  and  Mirie, 6  and  the  numbers  in  square 
parentheses  are  interpolations  of  the  numerical  results  of  Fenton 
and  Rienecker.7 

Columns  6  and  7  indicate  how  far  away  from  the  center  of  the 
collision  the  minimum  amplitude  occurs  (see  Fiq.  3),  and  by  hew 
much  this  minimum  amplitude  is  less  than  the  initial  amplitude. 
Fenton  and  Rienecker  report  their  numbers  as  a  final  chanqe  of 
amplitude,  but  it  appears  that  their  computational  domain  does 
not  exceed  Ax^n  by  very  much  (see  their  Fig.  4).  Entries  to 
columns  8  and  9  come  from  analyzing  results  at  very  late  times. 
The  phase  is  the  time-lag  of  an  outgoing  wave  extrapolated  from 
late  times  to  the  collision  time.  The  range  of  amplitude  changes 
is  set  down  to  account  for  the  failure  of  the  numerical 
calculations  to  properly  conserve  energy.  The  smaller  number  is 
the  computed  amplitude;  the  larger  number  takes  away  from  the 
solitary  waves  all  of  the  energy  added  by  the  numerics,  using 
energy-amplitude  curves  of  Longuet-Higgins  and  Fenton. 11 


«  3  4  1 


DIMENSIONLESS  DISTANCE 

Figure  1  —  Head-on  collision  of  an  identical  pair  of  solitary  waves  of 
amplitude  0.386.  Computations  are  carried  out  in  a 
computational  channel  of  3200  ecrually-spaced  grid  points  using 
Ax  *  At  »  1/8.  The  profiles  shown  are  eoually-spaced  in  time, 
starting  at  the  bottom,  with  the  vertical  spacing  between 
profiles  set  so  that  disturbances  that  propagate  at 


'dimensional 


'gH  would  be  aliened  along  ±  45  in  the  fioure. 
o 


Figure  2 


MAXIMUM  AMPLITUDE 


DIMENSIONLESS  TIME 

Figure  3  —  Time  history  of  a  crest,  (a)  Amplitude:  The  crest  amplitude 
goes  far  off  scale  during  the  collision  (to  0.893).  The  small 
wiggles  in  the  vicinity  of  time  ■  150  are  not  an  artifact  of 
the  machine >dravn  graphics,  but  may  be  a  numerical  artifact; 
other  second  order  versions  do  not  generate  them.  All  versions 
asymptote  to  amplitudes  less  than  initial  by  0.0009. 
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COLLIDING  WAVE  PHASE 


Tima  history  of  «  exoat.  (b)  Phase:  The  data  are  analyzed  to 
give  the  speed  and  a  reference  location  of  an  Incoming  solitary 
wave.  A  tracer  soring  with  this  speed  and  lined  up  at  the 
reference  location  la  defined.  The  plot  shoes  the  difference 
betveen  the  position  of  a  wave  creat  and  the  position  of  the 
tracer.  Between  the  tiaws  of  earleue  and  minimum  on  the  curve 
the  create  are  merged  at  x  ■  200. 
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